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Models with three active neutrinos and one sterile neutrino can naturally account 
for maximal oscillations of atmospheric neutrinos, explain the solar neutrino deficit, 
and accommodate the results of the LSND experiment. The models predict either 
i/ e —* v T or Of —- u s oscillations in long-baseline experiments with the atmospheric 
5m 2 scale and amplitude determined by the LSND oscillations. 


1 Introduction 

When neutrino flavor eigenstates Vf are not the same as the mass eigenstates 
Vi, e.g., for two neutrinos, 

Vf = cos Ov i + sin Ov 2 , Vf = — sin 0v\ + cos dv 2 , (1) 

then neutrinos oscillate. The vacuum oscillation probabilities are 


P{vf Vf>,L) = A sin 
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“appearance”, ( 2 ) 

“survival” , (3) 


where A = sin 2 20, 5m 2 = m\ — m 2 ; L is the path length and E is the neutrino 
energy. 

There is mounting experimental evidence for neutrino oscillationsl from 
the atmospheric neutrino anomaly (vacuum oscillations), the solar neutrino 
deficit (matter or vacuum oscillations), and the LSND experiment (vacuum 
oscillations). Each can be explained by oscillations of two flavors. However, 
three independent Sm 2 are required, but there are only two independent 5m, 2 
from v e , Vy,, and v T . If all observed oscillation effects_are real, a way out is 
oscillations to both active and sterile neutrino flavorstm 

Sterile neutrinos have no electroweak interactions (e.g., Z 7 4 v s v s ) and 
thus evade accelerator constraints. However, active <-> sterile oscillation in 
the early universe would lead to N„ = 4 neutrino species by the time of big 
bang nucleosynthesis (BBN), which is inconsistent with an N v % 2.6 bound 
based on low deuterium abundance but allowed by conservative estimates that 
N v < 4.5. If there exists a lepton number asymmetry L„ = (n v — n s )/n 7 > 
7.5 x 10 -5 at the epoch with temperature T < 10-20 MeV, then the appearance 
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of sterile neutrinos in the early universe can be suppressedu In the following, 
both tightly constrained (<5m 2 s A/ s < 10~ 7 eV 2 ) and unconstrained (Af s ~ 1) 
active <-> sterile oscillation possibilities are considered. 

2 The Data 

LSND The Los Alamos experiment studied z/ M —> v e oscillations from of 
decay at rest and —> i/ e from of 7r + decay in flight. The results, including 

restrictions from BNL, KARMEN and Bugey experiments, suggest —> i/ e 

oscillations with 

0.3 eV 2 < <5to 2 snd < 2.0 eV 2 , A LS nd « 4 x 10 " 2 to 3 x 10” 3 . (4) 

Atmospheric Cosmic ray interactions with the atmosphere produce 7 r-mesons 
and the decays 7 r — > jjlv and // —> v e ev^ give and v e fluxes in the approximate 
ratio (v,j, + P M )/(i/ e + v e ) ~ 2 . Measurements of R = (V M /iV e ) da t a / {N fl /N e ) M c 
for ~ 1 GeV find values of i? ~ 0.6. In the water Cherenkov experiments 
the single rings due to muons are fairly clean and sharp, while those from 
electrons are iuzzy to do electromagnetic showers. The Super-Kamiokande 
measurementof R versus the zenith angle 9 axe shown in Fig. la for sub-GeV 
and multi-GeV energies. As suggested long agoQ the data are well described by 
Vy, —> v T or —> v s oscillations with 5m\ TM w 5 x 10 _ 3 eV 2 and Aatm ~ 1, 
as shown by the dotted histograms in Fig. la. The relation of the path length 
L to the zenith angle 6 is displayed in Fig. lb. For sub-GeV neutrino energies, 
L/E is large at cos 9 < 0 and the oscillations average, P{v^ —> v^) « 0.5. At 
multi-GeV energies, L/E is large at cos# = —1 and P{v^ —> v^) « 0.5; also 
L/E is small at cos 9 = +1 and P{y^ —> v^) ss 1. The separate distributions of 
/i-like and e-like events versus the zenith angle establish that the anomalous R- 
ratio is due to a deficit of upward //-like events. The allowed ranges of — > v T 
oscillation parameters are summarized in Fig. 2. 

Solar Three types of solar v e experiments, (i) v & capture in Cl [Homestake], 
(ii) v e e —> v e e [Kamiokande and SuperKamiokande], (iii) v e capture in Ga, 
measure rates below standard model expectations. The different experiments 
are sensitive to different ranges of solar E v . There are three regions, of oscilla¬ 
tion parameter space that can accommodate all these observations 0 
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Figure 3 illustrates these parameter regions for the solar solutions along with 
the regions for the atmospheric and LSND oscillation interpretations. The 
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Figure 1: (a) The zenith angle dependence of R for sub-GeV and multi-GeV atmospheric 
neutrino samples from Super-Kamiokande; from Ref. pi (b) Path length versus zenith angle. 


solar v oscillation solutions will eventually be distinguished by use of all the 
following measurements: (i) time-averaged total flux, (ii) day-night dependence 
(earth-matter effects), (iii) energy spectra (electron energy in ve —> ve events), 
(iv) seasonal variation, and (v) the neutral-current to charged-current event 
ratio (SNO experiment). 
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Figure 2: Confidence intervals for sin 2 20 and Am 2 based on a % 2 fit to Super-Kamiokande 
atmospheric neutrino data. The solid line is 90% CL, the dashed line is 99% CL. From 
Ref. §! 



Figure 3: The three allowed two-neutrino solar solutions for v e —> v T oscillations. The 
corresponding regions for z/ e —>► Vs oscillations are similar to the v e —*► Vt case. 
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3 4-Neutrino Models 


Table 1 shows the options for oscillation solutions to all data. The preferred 
mass spectrum is two nearly degenerate mass pairs separated by the LSND 
scale, as displayed iix. Fig. 4. Here too > mi for the case of matter oscillations 
so that v e is resonant] in the Sun. The alternative of a 1+3 mass hierarchy with 
one heavier mass scale separated from three lighter, nearly degenerate states 
is disfavored when the null resultS|-|Of reactor and accelerator disappearance 
experiments are taken into accounts 

SOLAR 


LSND ATM SAM LAM VLW 



V e A 

Ve -» V T □ □ □ 


Table 1: Four-neutrino oscillation possibilities. 
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Figure 4: Neutrino mass spectrum showing which mass splittings are responsible for the 
LSND, atmospheric, and solar oscillations. 


Consider first the case with v, 
A neutrino mass matrix of the for: 


for the atmospheric oscillations. 
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can reproduce the three observed Sm 2 , the amplitudes Alsnd and Asol, and 
naturally give A atm = 1. Here the e^ are all small compared to unity. The 
values of ei and £2 determine which of the three solar solutions is realized. The 
mixing matrix corresponding to this mass matrix is 


V s ^ 


( cos 9 — sin 0 0 0 ^ 


( Vo ^ 

V & 


sin 0 cos 9 -4=e 3 ^e 3 
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The vacuum probabilities in this model are 


LSND 

ATM 

SOLAR 


P{y» <-► Ve) = 4 4 sin 2 

tj, \ ■ 2 (m 2 e 4 L\ 
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P(v e <->• v a ) ~ 2 2 2 sin 2 Asol 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


where A S ol = 5m\ OIj L/E, <5to 2 ol ~ 4e^/ef, 4e|/(ef + el), and 1, for SAM, 
LAM, and VLW, respectively. Matter effects must be included in the solar 
SAM and LAM solutions. 

By construction this model has effective two-neutrino oscillation solu¬ 
tions for LSND, ATM, and solar phenomena. The model makes a number 
of predictionsB 

(i) Neutrinoless double-/? decay vanishes at tree level because M VeVe = 0. 

(ii) The neutrino mass spectrum is 7713,7714 ~ 1.4 eV, mo ~ 2 x 10 -3 eV, 
TOi ~ 4 x 10~ 6 eV. There will be no measurable effect at the endpoint of tritium 
beta decay if v e is primarily associated with the lighter pair. 

(iii) The hot dark matter contribution to the mass density of the Universe 
is = ^2(m v /93 eV)h~ 2 , where h ~ 0.65. The SLOAN Digital Sky Survey is 
expected to.have sensitivity down to m v = 0.2-0.9 eV for two nearly degenerate 
neutrinosPdl which covers the interesting range from LSND. 

(iv) In the SNO solar experiment, both CC and NC event rates would be 
suppressed, with NC/CC = 1, if v e —7 v s is the solar solution. 

(v) In reactor experiments the v e disappearance P(v e —► v e ) = 1 — 
Alsnd sin 2 Alsnd with Alsnd ~ 2.5 x 10~ 3 is not detectable. For example, 
the CHOOZ experiment sensitivity is A > 0.2 for 6m 2 <1 10 _3 eV 2 . 
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(vi) Long-baseline experiments with L/E « 10-10 2 km/GeV could mea¬ 
sure P(v^ «-> v T ) — sin 2 Aatm and confirm the atmospheric oscillation result. 
In addition, the prediction of new oscillations 

P(v^<r^v e ) ~ a(2 - sin 2 Aatm) (10) 

P(y e ^v T ) ~ a sin 2 Aatm (11) 

with a = Alsnd/4 « 10 -2 to 10 -3 could be tested. For this purpose intense 
neutrino beams are required. The MINOS experiment (Fermilab to Soudan) 
could confirm the v f , —■> v T oscillations and test the —+ v e prediction, pro¬ 
vided that 5 toatm > 2 x 10 -3 eV 2 . 

In the future, a special purpose muon storage ring could provide high 
intensity neutrino beams with welLdetermined fluxes that could be directed 
towards any detector on the earthlB It could be possible to store ~ 10 21 /i + or 
fi~ per year and obtain ~ 10 20 neutrinos from the muon decays. Oscillations 
give “wrong sign” leptons from those produced by the beam. For example, 
H~ decays give v e and fluxes so detection of /x + , e~ , r ± leptons tests for 
Ve —> Vy,{v T ) and v /t — ¥ v e {y T ) oscillations. Taus can be detected via their 
r —> [i decays and the r-charges so determined to distinguish v^ —> v T and 
v e —> v T oscillations. The ranges of oscillation parameters that could be tested 
in such long-baseline experiments is illustrated in Fig. 5. 

4 Alternative 4-Neutrino Mixings 

The preceding model assumed large v lt <-> v T mixing as the explanation of 
the atmospheric data. The alternative scenario with large v^ v s mixing is 
obtained in a straightforward manner by interchange of v s and v T labels. Then 
the long-baseline predictions are v^ disappearance 

P{vf, v s ) ~ sin 2 Aatm (12) 

and v /J: <-> v e appearance 

P{y^ ++ v e ) ~ a(2 - sin 2 A AT m) , (13) 

A more general scenario could have large v,, mixing with a linear combination 
of v s and v T . 

5 Summary 

A simple mass matrix for four neutrinos with strategically placed zeros can 
accommodate the indications for neutrino oscillations from the LSND, atmo- 
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Figure 5: Predicted region in the effective 8m 2 - sin 2 20 parameter space for v e —> v T oscil¬ 
lations in the four-neutrino model (solid rectangle), which is determined by of the LSND 
—> v e oscillation amplitude and the atmospheric neutrino v^ > v T oscillation 8m 2 scale. 

The dotted curves show the potential limits on —> v e ,v T oscillations from the MINOS 

experiment and the dashed curves show the potential limits on i/ e , —► i>v oscillations that 

can be set by neutrino beams from an intense muon source at Fermilab to detectors at the 
SOUDAN and GRAN SASSO sites for muons with energy of 20 GeV. Also shown are the 
parameters for the solar v e —* v s small-angle MSW oscillation. 


spheric, and solar data. The two principle variants of the model have oscilla¬ 
tions of 


—► v T or v s ATM (14) 

v e —► v s or v T SOLAR (15) 

with Vp —> v e for LSND. The predictions for long-baseline experiments are 
oscillations with the scale with maximal amplitude in the channel 

Vp —► v T or v s and amplitudes ~ 10 -2 to 10 -3 in the channels v e —> v T or v s 
and i/fj, —> u e . 
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